Quantum Hall states of bosons in rotating anharmonic traps 



Nicolas Rougerie 

Universite Grenoble 1 & CNRS, LPMMC (UMR 5493), B.P. 166, 38 042 Grenoble, France. 



CO 

o 

CN 

X> 
<D 
IX, 

CO 



wo; 

a* 
3 



Sylvia Serfaty 

UPMC Univ Paris 06, UMR 7598 Laboratoire Jacques-Louis Lions, Paris, 
F-75005 France ; CNRS, UMR 1598 LJLL, Pans, F-75005 France and 
Courant Institute, New York University, 251 Mercer st, NY NY 10012, USA 

Jakob Yngvason 
Faculty of Physics, University of Vienna, Boltzmanngasse 5 and 
Erwin Schrddinger Institute for Mathematical Physics, Boltzmanngasse 9, A-1090 Vienna, Austria. 

(Dated: January/30/2013) 

We study a model of bosons in the lowest Landau level in a rotating trap where the confinement 
potential is a sum of a quadratic and a quartic term. The quartic term improves the stability of 
the system against centrifugal deconfinement and allows to consider rotation frequencies beyond 
the frequency of the quadratic part. The interactions between particles are modeled by a Dirac 
delta potential. We derive rigorously conditions for ground states of the system to be strongly 
correlated in the sense that they are confined to the kernel of the interaction operator, and thus 
contain the correlations of the Bose-Laughlin state. Rigorous angular momentum estimates and trial 
state arguments indicate a transition from a pure Laughlin state to a state containing in addition a 
giant vortex at the center of the trap (Laughlin quasi-hole). There are also indications of a second 
transition where the density changes from a flat profile in a disc or an annulus to a radial Gaussian 
confined to a thin annulus 

PACS numbers: 67.85-d,05.30. Jp,03.75.Hh 
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I. INTRODUCTION 

One of the most striking phenomena in condensed 
matter physics is the Fractional Quantum Hall Effect 
(FQHE) for charged fermions in strong magnetic fields 
that still, after decades of research, poses many challeng- 
ing questions [l|-|3|. 

The formal similarity between the Hamiltonian of a 
rotating Bose gas and that of a 2D electron gas in a 
magnetic field suggests the possibility of studying bosonic 
analogues of the FQHE in cold quantum gases set in rapid 
rotation 0-11]. In particular, strongly correlated phases 
analogous to the fermionic Laughlin state 0, [l(| are pre- 
dicted to occur in rotating boson clouds. More specifi- 
cally, a phase transition from a Bose-Einstein condensate 
(BEC) with a vortex lattice to a strongly correlated state 
should happen at total angular momentum oc N 2 (filling 
factor of order 1), where JV > 1 denotes the particle 
number. This phase transition has not been observed 
yet because it is extremely difficult to reach such high 
angular momenta in the laboratory. 

The reasons for this difficulty can be understood as fol- 
lows. The usual experimental setup consists of bosonic 
atoms in a rotating trap modeled by a harmonic poten- 
tial. Once the usual approximations have been made, the 
relevant Hamiltonian boils down to 



H = J2 



■g^SQci-Xj) (LI) 



3=1 



i<j 



operating on wave functions ^(x\ 



,x N ), x, e 



the lowest Landau level (LLL) of a two-dimensional mag- 
netic Hamiltonian. Here u> > is half the difference be- 
tween the squares of the frequency of the harmonic trap 
and of the rotation velocity, measured in units of the 
latter, and g > is the coupling constant of the inter- 
action. To obtain correlated states, one wants to favor 
the interaction, i.e., consider small values of oj / g. Clearly 
this corresponds to a singular limit where the gas is no 
longer confined against centrifugal forces. This is the 
main source of difficulty that has to be overcome to create 
the Bose-Laughlin state (see [Tljj for a more quantitative 
discussion and |12l . |13| for state of the art experiments in 
this direction). 

A natural proposal [1, [HI, E3 is to add a confining term 
that should be stronger than quadratic (i.e., stronger 
than the centrifugal force) . In the present paper we inves- 
tigate asymptotic properties of ground states after such 
an addition to the Hamiltonian for large particle num- 
bers and when the parameters are tuned to favor strong 
correlations. A simple model for the additional term is 
k\xj\ 4 with k > 0, that has already been considered in 
the context of rotating BECs [HI, [lj| . Besides the quest 
for the Laughlin state, our motivation is also to explore 
new phenomena arising from this modification of (|I.ip . 

A new problem that has to be tackled is the lack of 
commutativity of the anharmonic potential term, pro- 
jected onto the lowest Landau level, and the interaction. 
This implies that the Laughlin state is not an exact eigen- 
state of the modified Hamiltonian as it is for the Hamil- 
tonian in a purely harmonic trap. Nevertheless we prove 
that Laughlin-like correlations can be achieved asymp- 
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totically by a suitable choice of the parameters ui and k 
in dependence of N. Moreover, we argue that the ground 
state exhibits a new strongly correlated phase if uj is neg- 
ative and the ratio \u>\/k of order N or larger. In this 
phase the density is depleted in a hole around the center 
of the trap while it has the features of an incompress- 
ible quantum fluid in an annulus around the hole. In 
fact, there are clear indications of a further transition at 
MA ^ N 2 with the density concentrated in a thin an- 
nulus of large radius with a Gaussian profile in the radial 
variable, while for \u)\/k significantly smaller that iV 2 the 
profile in the annulus is essentially flat. 

Our analysis relies on the one hand on energy esti- 
mates, making use of a representation of the anharmonic 
term through angular momentum operators, and on the 
other hand on quantitative estimates for the one-particle 
density of strongly correlated trial states. These states 
have the form of a Laughlin state modified by a factor 
corresponding to a vortex of high angular momentum 
at the origin ('Laughlin times giant vortex' or Laugh- 
lin quasi- hole). Mathematically the energy estimates are 
much simpler than the second part of the analysis and al- 
ready lead to criteria for strong correlations in an asymp- 
totic limit and estimates on the angular momentum of 
the ground state. The estimates on the density, on the 
other hand, are crucial for a physical interpretation of the 
two phase transitions indicated by the energy estimates. 
Moreover, the density estimates imply an improvement 
of the energy estimates. 

The intuition behind the density estimates is Laugh- 
lin's plasma analogy 0, [l(|: The trial states are seen 
as defining Gibbs measures of a 2D classical Coulomb 
gas. The one-body densities appear as combined mean 
field/low temperature limits for the Coulomb gas of the 
type studied e.g. in [l7l - [20l | using compactness arguments. 
For quantitative error estimates, a new method for the 
study of the mean-field limit is needed. In the present 
paper we outline this method but refer to a companion 
paper [2l[ for full details. 

The depletion of the density at the center of the trap 
and absence of vortices in the bulk for appropriate val- 
ues of the parameters is a feature common to the so- 
called 'giant vortex' phase of uncorrelated states in the 
Gross-Pitaevskii (GP) approximation. The transition 
from a vortex lattice to a giant vortex phase in anhar- 
monic traps has been studied by various methods over 
the past 10 years, e.g., in |22h27| . and was recently es- 
tablished with full mathematical rigor (see 15[ for refer- 
ences). We stress, however, that this transition is phys- 
ically and mathematically quite different from the one 
studied here. 

The rest of the paper is organized as follows: In Section 
IIVI below we state our results on the energy estimates 
and their consequences in the form of three theorems. 
The energy upper and lower bounds match to leading 
order, and the change of optimal trial functions when 
the parameters are varied indicates a transition from one 
type of strongly correlated ground state (Laughlin type) 



to another (Laughlin times giant vortex). The precise 
statements of the results requires some preliminary dis- 
cussions in Section Ql] of the properties of the Hamilto- 
nian (jl.ip regarded as an operator on a Bargmann space 
of analytic functions. In Section [IIII we introduce the an- 
harmonic addition to the potential as well as its alternate 
definition through the square of the one-particle angular 
momentum operator in the lowest Landau level. The the- 
orems stated in Section IIVI are proved in Section [V] In 
Section IVT1 we discuss the determination of the one-body 
density in the strongly correlated trial states considered 
and its implications for the interpretation and improve- 
ment of the results of Section IIVI A final Section IVIII 
summarizes our findings and discusses their relation to 
the recent paper [TTj . 

More details on our approach can be found in the com- 
panion paper (2lj . 



II. THE STANDARD HAMILTONIAN 

It is well known [H, |2!| (see also [3(J HH for more 
references) that the Hamiltonian (|I.1[) in the LLL can be 
regarded as an operator on the Bargmann space Sat of 
symmetric analytic functions F of zt, ■ . . , Z/v with £ C 
such that 

/ \F(zi, . . . ,zjv)| 2 exp (-Y,7=i\ z j\ 2 ) dz i ' ' ' dz N < oo. 

Jc N 

(HI) 

Here dz denotes the Lebesgue measure on M 2 that is iden- 
tified with the complex plane C in the usual way. On Bn 
the angular momentum operator of the i-th particle is 
Li — Zid Zi and the contact interaction S(zi — Zj) is given 
by the bounded operator 



8 l3 F{. 



■) 



= ^ ...,!(* + ... , !(% + Zj), ...) {11.2) 

as noted first in |29j. Indeed, using the analyticity of F, 
it is easy to sec that 

(F )( %F)= f |F(...,z,...,z,...)| 2 exp(-2|z| 2 )dz 

Jc N 

xexp(-E^ J l^| 2 )n fe ^, i ^ ! (H.3) 

where (•) is the scalar product on corresponding to 
(III. II) . Using that zd z ■<-» \z\ 2 — 1 in this scalar product 
the operator (|I.1[) can, apart from a constant additive 
term Nuj, be written as 



H N = uj£ n + gl 



N 



with 



N 



(II.4) 



(II.5) 
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We denote by Hn the space of wave functions of the 
form %(zi,...,z N ) = F(zi, . . . , z N ) exp(- £V =1 M 2 /2) 
with F £ Bn- This is a subspace of L 2 (M. 2N ) and clearly 
isomorphic to Bn- Note, however, that the angular mo- 
mentum operator Li on Hn is Zid Zi — z~id Zi rather than 
Zid Zi as on Bn- To distinguish these two realizations of 
states in the LLL we shall denote functions in Bn by 
roman letters and in Hn by greek letters. 

We remark that the reduced Hamiltonian (|II.4[) is ob- 
tained by restricting states to the LLL associated with 
the rotational frequency (that we fix equal to 1), which 
leads to p.f|) . and then using the substitution zd z 
\z\ 2 — 1. Another possibility (see [1, HO, UK and refer- 
ences therein) is to restrict states to the LLL associated 
with the harmonic trap frequency, which directly leads 
to the model (|II.4[) but with w equal to the difference be- 
tween the two frequencies (rather than half the difference 
between the squares) . Since in the regime of interest the 
two frequencies are close, the two approaches are equiva- 
lent. Ours has the advantage of better emphasizing what 
the residual effective potential is, which is natural when 
we consider the additional anharmonic term below. 

An essential feature of the Hamiltonian (|II.4[) is that 
the operators Cn and In commute. The lower bound- 
ary of (the convex hull of) their joint spectrum in a plot 
with angular momentum as the horizontal axis is called 
the yrast curve (see, e.g., Q and [3(| for plots show- 
ing its qualitative features). As a function of the eigen- 
values L of Cn the yrast curve I(L) is decreasing from 
7(0) = (47r)- 1 7V(7V - 1) to I(N(N - 1)) = 0. The mono- 
tonicity follows from the observation that if a simulta- 
neous eigenfunction of Cn and In is multiplied by the 
center of mass, (z% + • • • + zn)/N ', the interaction is un- 
changed while the angular momentum increases by one 
unit. 

For a given ratio oo/g the ground state of (|II.4|) (in 
general not unique) is determined by the point (s) on the 
yrast curve where a supporting line has slope —oj/g. 



A. Fully correlated states 



For L < N the 



ground state of (|II.4[) is explicitly 
known (see [32j or |8| and references cited there) while 
for large TV and L <C TV 2 a Gross-Pitaevskii description 
with an uncorrelated ground state is asymptotically cor- 
rect [3lj]. Results about the Gross-Pitaevskii functional 
with states restricted to the lowest Landau level may be 
found in [33T - l36j | . 

For L = TV (TV — 1) the unique ground state of In with 
eigenvalue is the bosonic Laughlin state whose wave 
function in Bn is the symmetric polynomial 



Fi*m(zi, • • • , Zn) = const. \\{z l - ZjY 



(II, 



i<3 



or, equivalently, uyHn 

^Lau^i,. . . ,z N ) = const. U^i-z,-) 2 exp(-X)fcl 2: fc| 2 / 2 )- 

(II.7) 

More generally, the null space Kerljv of the interaction 
operator consists of functions of the form 

F(zi, z N ) = G(zi, ZAr)F Lau (zi, ...,z N ) (11.8) 

where G is a symmetric analytic function such that F is 
square integrable w.r.t. Y\ k exp(— |zfc| 2 )dzfc. We shall call 
states of this form fully correlated states. Their angular 
momentum spectrum is contained in L > TV(TV — 1). 

A common interpretation of the form (|II.6|) is that each 
particle "binds" a vortex a degree 2 (see e.g. 0, H and 
references therein). Indeed, since LLL wave functions are 
analytic, putting a zero at each particle in order to cancel 
the interaction energy requires a non trivial quantized 
phase. 



B. Spectral gaps 

States with a given eigenvalue L of the total angular 
momentum correspond to symmetric, homogeneous poly- 
nomials of degree L. This is a finite dimensional space 
and for all L it contains states with strictly positive inter- 
action energy. (Take for example F — (const.) J2k z k-) 
Hence the spectral gap 

A(L) = smallest eigenvalue > of In \{c n =l\ (H.9) 

is strictly positive for all L. Note that A(L) depends on 
TV besides L. 

The trial state argument proving the monotonic de- 
crease of the yrast curve also applies to A(L), and proves 
that it is decreasing with L. Numerical studies [El. l37l - [39l | 
suggest that A = inf L A(L) = A(N(N - 1) > 
and it is also believed (see the discussion in [30]) that A 
stays of order 1 and bounded away from zero as JV — > oo. 
If this spectral gap conjecture is true, the statements of 
Theorems [5] and [3] below can be simplified, but to our 
knowledge the conjecture has not yet been proved. We 
shall therefore not rely on any such assumptions but state 
our results in the sequel in terms of A(L). 

According to numerical diagonalizations d, H3, IH[ , it 
seems that the yrast state corresponding to the interac- 
tion energy A(N(N— 1) — 1) is a center of mass excitation 
of the the state corresponding to A(N(N— 1) — TV), which 
would imply A(TV(TV - 1) - 1) = A(TV(TV - 1) - TV). The 
latter quantity is associated with the energy cost to cre- 
ate a Laughlin quasi-particle 0, in the system. We 
do not know a proof that this picture is correct, however. 



III. ADDING AN ANHARMONIC POTENTIAL 

The limit u> —¥ keeping u > is experimentally very 
delicate. For stability, but also to study new effects, we 
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consider now a modification of the Hamiltonian where 
one adds a quartic term to the potential: 



N 



(III.1) 



with a new parameter k > 0. The expectation value of 
the energy for a normalized ^> — F ■ exp(— J^. \zi\ 2 /2) £ 
T-Ln is given by the functional 

£[*] = y V^tizjp^z) + g{F,l N F) (111.2) 

where 

VL, fc (z) =w|z| 2 + fc|z| 4 (III.3) 
and pq, is the one-body density with the normalization 



p\&(z) dz = N. 



(III.4) 



In particular, the energy is bounded below by N min V Ut k 
that is finite even if id < 0, provided k > 0. For the fully 
correlated states (|II.8|) the second term in (|III.2|) van- 
ishes and the expectation value of the energy is entirely 
determined by the one-particle density. 

Using the correspondence zd z ■<-» \z\ 2 — 1 and (zd z ) 2 O 
\z\ A — 3|z| 2 -|-l the modified Hamiltonian on the Bargmann 
space Bn can be written (up to the additive constant 
N(uj + 2k) that will be dropped in the following) as 



N 



H' N 



(u + 3k)C N + k^L 2 + gX N . (III.5) 



1=1 



For the original Hamiltonian (|II.4[) . without the anhar- 
monic addition to the potential, the Laughlin state is 
an exact eigenstate with L = N(N — 1) and energy 
uiN(N— 1). For k ^ 0, however, the Laughlin state is not 
an eigenstate of H' N because Yli=i does not commute 
with In . On the other hand Cn still commutes with H' N 
and we may consider the ground state energy Eq(L) of 
H' N restricted to the subspace where Cn = L. The unre- 
stricted ground state energy min/, Eq(L) will be denoted 
by Eq, a ground state (that might not be unique) by ^q 
and its angular momentum by Lq. 



IV. ASYMPTOTIC PROPERTIES OF GROUND 
STATES 

We now state our main results concerning the angu- 
lar momentum Lq, the confinement to the space of fully 
correlated states, and the energy of a ground state ^q 
of pil.5[) in the limit when N — » oo and at the same 
time uj, k — > 0. Note that "Jo can be assumed to have a 
definite angular momentum because Cn commutes with 

(Him . 



Theorem 1 (Angular momentum) 

In the limit N — > oo, w, k — > the angular momentum 
Lq of a ground state of H' N satisfies 

• Ifid> -2kN 



L Q < 2N 2 
Ifid< -2kN but \uj\/k < N 2 

\LQ-L qh \ < v/3iV 2 (l+o(l)), 

where 

\id\N 



Lqh '■- 



2k 



Ifcu< -2kN and \id\/k > N 2 

|£o-iqh| < V3L^ 2 7V(l + (l)). 



(IV.l) 
(IV.2) 
(IV.3) 

(IV.4) 



In particular, Lq/L^ — > 1 in both cases (|IV.2[) and 
(ITOll . 

The next theorem provides criteria for a ground state 
to be asymptotically fully correlated, i.e., in the kernel of 
the interaction operator. In the statement the following 
gaps in the spectrum enter: 

A! = A(2iV 2 ), A 3 = A(L qh + V3(N 2 )), 

A 4 = A(L qh + y30L^ 2 AO). (IV.5) 

If the spectral gap conjecture holds, all these gaps can 
be replaced by the minimal gap A > independent of N 
and the parameters. 

Theorem 2 (Criteria for full correlation) 

Let PkciIXn) 1 - denote the projector on the orthogonal 
complement of the kernel of In . We have 



KKcr(I JV ) i *o|| i 2( R2 ]V) -> 



(IV.6) 

if one of the following 



in the limit N — > oo, id, k 
conditions holds: 

Case 1. id > and (g Ai)" 1 • (idN 2 + kN 3 ) ->■ 0. In this 
case, for N large enough, 

\\PK Cli i N )^o\\ 2 < (.9 A^- 1 ■ <3kN 3 (l + o(l)). (IV.7) 

Case 2. -2kN < u j < and (g A 1 )- 1 -(N(id 2 /k)+cdN 2 + 
kN 3 ) — s- 0. Again (|IV.7|) holds under this condition for 
N large. 

Case 3. id < -2k N , \id\/k < N 2 and (g A3)- 1 • kN 3 -> 
0. Here 



\\P Kci iT N )^o\\ 2 < (g ^-'-(Sk N 3 + -\id\N ) (l+o(l)). 

(IV.8) 

Case 4. id < -2k N , \w\/k > A^ 2 and (g A4)- 1 ■ \id\N -> 
0. Here 

II^K CT( i N )^o|| 2 < (9 A4)- 1 • l\td\N (l + o(l)). (IV.9) 
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Both Theorems 1 and 2 are consequences of the energy 
bounds stated in the following Theorem 3. The "cases" 
referred to are the same as in Theorem 2. 

Theorem 3 (Energy bounds) 

The ground state energy Eq satisfies the following bounds: 
Cases 1 and 2 

{ujN 2 + kN 3 )(l - o(l)) <E Q < ujN 2 + AkN 3 . (IV.10) 

Case 3 

-^(1+ (1)) < E < -^L+(skN 3 + 1\ U \N) (1+o( 

(IV.ll) 

Case 4 

u 2 N 



m > are often referred to as 'Laughlin quasi-holes'. 
They already appeared in Laughlin's seminal paper Q. 
Note that we will take a large value of m in the sequel, 
in which case the denomination 'Laughlin times giant 
vortex' seems more appropriate. 

Since the trial states (|V.2j) all have a zero interaction 
energy, we only have to estimate their potential energy. 
We have 



^^F^ = (£nG)Fl &u + G(£/vi<Lau) 



(Nm + N(N- 1))F£K (V.3) 



1)) 



4fc 



(l + o(l)) < E Q < 



uj 2 N 
4fc 



Mjv(i + (i)). 

(IV.12) 



Note that the present paper is only concerned with 
properties of the ground states of the system. The usual 
(heuristic and numerical) arguments supporting the pic- 
ture of fractionally charged excitations of the ground 
states [H, IH carry over with little modifications in our 
setting. They are not on the same level of mathematical 
rigor as the above theorems however, and we have no rig- 
orous arguments allowing to identify the excited states of 
our model. 



PROOFS OF THEOREMS 1-3 



The proofs of Theorems 1-3 rely on two ingredients: 

• A lower bound for the energy at fixed angular mo- 
mentum L: 



E (L) > (uj + 3k)L + k—. 



(v.i) 



• An upper bound for the energy of suitable trial 
functions. 



The proof of (|V.1|) is quite simple: Dropping the non- 
negative interaction term in (|III.5[) one has to minimize 

(F,((oj + 3k)C N + fcE l =i i ?)^) over F e B N with 
CnF = LF. The bound follows from the operator in- 
equality Yli^i > jriY^i^i) 2 that holds because L, and 
Lj commute for any i,j. 

The upper bounds on Eq can be derived as follows. 
As trial functions in the Bargmann space Bn we take 



qh 



JV 



G{z x , . . .,z N ) = c N<m JJ; 



(V.2) 



with CN, m such that F^' is normalized. One recovers 
the pure Laughlin state for m = and the states with 



Moreover, 



N 



N 



i=l 



Nm 2 + 2mN(N - 1) + N(GL 1 F LeiU , GLiF Lau ) (V.4) 

where the symmetry of -FLau has been used for the last 
term which can be written as 



N 



J |LiF Lal J 2 |G| 2 exp(-^N 2 );Qd^ 



(V.5) 



We now note that the monomials 2™ are or- 
thogonal for different n's in the L 2 -scalar product 
(. , .)q on analytic functions defined by the measure 
|G| 2 exp(— J^i \ z i\ 2 ) Hi dzi because |G| 2 depends only on 
the absolute values of the z,-'s. Moreover, Liz™ = nz", 
and Fl &u is a polynomial of degree 2(N— 1) in z\. Hence 

N (GLiFlhu, GLiF Lau ) < AN{N - l) 2 (F Lau , F Lau ) G 
= 4A^-l) 2 (GF Lau ,GF Lau ) 
= 4iV(iV-l) 2 (V.6) 

where we has also used the normalization of our trial 
state. Inserting Eqs. (|V.3|) . (IV.4I) and (|V.6[) in the expec- 
tation value (-Fqb , H' N F^) leads to the upper bound 

(F^'KF^) < (ui + mmN + NiN-l)) 

+ k{Nm 2 + 2N{N -l)m + AN 3 ). (V.7) 

Optimizing over m gives the value 



J if u > -2kN 
iopt -<iM_ Ni{LJ< _ 2kN 



(V.8) 



Strictly speaking m opt must be an integer so we should 
take the integer value of the above. Since we are only con- 
cerned with orders of magnitudes this does not produce 
any change in the final result, and we ignore this detail. 
For m opt = this gives the upper bound in (jIV.10[) . 
For TO op t = (|o;|/2fe) — N > we obtain 



En < — 



LJ 2 N 

4fc 



+ 3fcjV J + -|w|iV (l + o(l)). (V.9) 
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The first term in (|V.9[) is equal to N times the mini- take values between and N(N — |) so we can conclude 



mum of the potential uj\z\ 



k\ 



The second term 



~ fcJV 3 dominates the third term ~ \uj\N ~ kNm op t for 
\iv\/2k < N 2 , while the converse holds for \u)\/2k > N 2 . 
Note also that L q ^ ps Nm op t for |w|/2fc ^> N. 

To proceed we note that by (jV.ip and the definition of 
the gap we have 

(oj + 3k)L Q + k^+gA(L )\\P Kcl{ x N ^\\ 2 <E Q . (V.10) 



For ui > —2kN it follows from (jV.101) and the upper 
bound in (|IV.10p that 



(cj + 3k)L a + k— < loN 2 + AkN 3 
i N - 



Cv.il) 

iV 

which implies 

L < 27V 2 (V.12) 

and hence 

A{L ) > A(2N 2 ). (V.13) 

For u) > 0, /(£) = (w + 3fc)i + kL 2 /N is nonnegative. 
Hence we see that in this case ||i : kor(ijv) i *o|| 2 — > if 

{gA(2N 2 )}- 1 {kN 3 +ujN 2 ) -> 0. (V.14) 

On the other hand, if ||-PKor(i N )- L ^'o|| 2 —> 0, then it is 
clear that Lq must be > N(N — 1) for AT large enough, 
and hence f (L ) > {ujN 2 + kN 3 ){l - o(l)). 

Altogether we have, for u) > and under the condition 
(|VT4l) . the bounds 

(wJV 2 + fcA^ 3 )(l - o(l)) < £ < wiV 2 + 4fcA^ 3 (l + o(l)) 

(V.15) 

as well as 



|P K cr(z N )-*o|| 2 <3{ 5 A(2Ar 2 )rW(l + o(l)). 



Consider next the case —2kN < u) < 0. As before, Eq. 
(jVXtl) implies (TVTT2|) and (fVU3l) . For the analogue of 
(IV. 141) we use that f(L) takes a minimal value at 

AT 

Lmin = ~2k^ + 3fc) = Lqh(1 + (V ' 17) 



with 

/(imin) 



.| (w + 3 fc) 2 ^-^(l + (l)). (V.18) 



Condition (|V.14[) for a fully correlated ground state thus 
gets replaced by 

{g A(2N 2 )}- 1 fkN 3 + ujN 2 + N^-) -> 0. (V.19) 



Moreover, (IV.19[) implies as before that Lq eventually 
becomes > N(N - 1). For -2kN < uj < 0, L min can 



that f (L ) > f(N(N- 1)) that leads again to (jVT5|) and 

(TvT6l) . 

For the case U) < —2kN we use 



f(L) = f(L miD ) + -(L-L miD ) 2 



(V.20) 



and obtain from (|V.10|) and (IV. L8|l . up to factors (1 
^(Lo~L min ) 2 + gA(L )\\P Kcl{lN) ^ \\ 2 < 3kN 3 . 



In particular, 



\L Q - L min \ < V3N 2 



(V.21) 



(V.22) 



and thus A(L ) > A(L min + V^V 2 )- T he con dition for 
full correlation, analogous to (|V.14I) and (|V.19[) becomes 



{gA(L min + VZN^j^kN 3 -+ 0, 



(V.23) 



indeed, 



||^Ker(Z„)^0|| 2 < \{g A(L min + V3N 2 )}^k N 3 . 

(V.24) 

Since /(L m i n ) = —NuP/Ak to leading order, the energy 
estimate is 

■ j: " N -{\ - o(l)) < E < + 3fcA^ 3 (l + o(l)). 



4fc 



4A- 



(V.25) 

The condition for fcA^ 3 <C uj 2 N/k is equivalent to |w|/fc 3> 
AT, i.e., m opt > AT. Note also that (|V.22|) implies 



\ L ~ L lh\ 1 



(V.26) 



if m op t > N. 

If —uj/k ^> A^ 2 as in Case 4, the error term 3kN 3 in 



(V.16) (|V.9[) is small compared to §|w|A r ~ kN m op t and we 



have 

E <-^- + ^\N(l + o(l)). (V.27) 
Thus (|V.2ip gets replaced by (as always, up to (1 + o(l)) 

^(Lo-L min ) 2 +gA(L )\\P Kel{lN) ^ \\ 2 < *\u\N. 



and (|VT22|) by 



(V.28) 



|£u - Aninl < V3 ( y ) A^ (V.29) 



1/2 



which implies pV.4j) and (|V.26|) for m op t > A^ 2 . 

The analogues of (j\T23| . dV^l) and (|V^5|> for m opt > 
A^ 2 are 

{<? A(L min + VSNL^y-^lN 0, (V.30) 
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IP, 



KcrtZjv^^O 



< 2^ A ^ 



and 



+ V3NL^ n )}~^\N, 
(V.31) 



w 2 ^ , xx uj 2 N 
-0—o(i))<Eo< 



4fc 



4/v 



| W |iV(l + (l)). 



(V.32) 



VI. THE ONE-PARTICLE DENSITY OF FULLY 
CORRELATED TRIAL STATES 

The main conclusions that can be drawn from the pre- 
vious sections may be summarized as follows: 

• The energy estimates are consistent with the pic- 
ture that the Laughlin state is an approximate 
ground state for positive u in Case 1, and also for 
negative u in Case 2, as long as \u)\/k < N. The 
angular momentum remains 0(N 2 ) in these cases. 

• When u) < and \w\/kN becomes large (Cases 
3 and 4) the angular momentum is approximately 
£ q h = 0(N\w\/k) > TV 2 , much larger than for the 
Laughlin state. 

• A transition between Cases 3 and 4 at \to\/k ~ TV 2 
is manifest through the change of the sublcading 
contribution to the energy of the trial functions. 
Its order of magnitude changes from 0(kN 3 ) to 
0(\lj\N) at the transition. 

We now want to substantiate this picture further by 
investigating the one-particle densities of our trial states, 
This will both clarify the physical meaning of the two 
transitions at \oj\/k ~ N and \oj\/k ~ TV 2 that are indi- 
cated by the energy considerations above, and also allow 
improvements of the energy estimates. The latter rely 
on the representation (|III.2|) of the energy and requires 
detailed estimates of the one-particle density. We note 
that the representation pil.2[) holds also for more general 
potentials than w|z| 2 + fc|z| 4 (e.g., anisotropic potentials) , 
where the approach via angular momentum presented in 
Section V would not apply. We shall here present the 
main results of the analysis of the density but refer to 
[2l| for details of the proofs. 

Already Laughlin's paper of 1983 j9j contains the im- 
portant idea that the determination of the probability 
density \^f\ 2 in his and related states for large iV is equiv- 
alent to a 2D electrostatic problem. We make use of this 
idea for our trial states whose wave functions in Hn are 

N 

..., ZN )= CN , m n zT n c* - e- ^ ^ |2/2 

(VI.1) 



where Cjv.m 1S a normalization constant. We denote 
(#!,..., zjv) by Z for short and consider the scaled prob- 
ability density 



fi N (Z) = N N *™{y/NZ) 



qh 



(VI.2) 



Then we can write 







Z JV,m ex P < 


Is- 


4^1og|^ 




i<j 




z n] m ex P 1 




_1 and 





(VI.3) 



H N (Z) = f:(\z j \ 2 - 2 -^lo g \z j \)-±J2^\ 

3=1 i<] 

The normalization factor (partition function) is 
Z N = J exp (-Lh n (Z^J d 2N Z. 



(VL4) 



In the sequel we discuss the properties of the density 
/ijv(^) in the iV — > oo limit. 



Mean field limit 



The Hamiltonian (|VI.4j) describes a classical 2D 
Coulomb gas in a uniform background of opposite charge 
and a point charge of strength 2m /N at the origin, cor- 
responding respectively to the \zj\ 2 and the — -Si log \zj\ 
terms. The probability measure (1n(Z) minimizes the 
free energy F(T, N,m) = — TlogZ^ for this Hamilto- 
nian at T = N' 1 . 

The N — > oo limit is in this interpretation a mean field 
limit at small temperature T — > 0. Thus, it is reasonable 
to expect that in the limit [in factorizes, /j.n ~ P® N > 
where the probability measure p on R 2 minimizes an 
appropriate mean field energy functional. In fact, us- 
ing compactness arguments, rigorous results of this type 
have been obtained in [18l - l2(j ] for related models. The 
computation of the energy via Equation (|III.2p . however, 
requires quantitative estimates of the errors in the ap- 
proximations that can only be obtained by a different 
method. A more elaborate analysis is developed in [Z(| 
to describe the next order fluctuations around p® N , how- 
ever it does not provide estimates in the form that we 
need here (cf. [2l|, Section 3] for more details). 

The mean field free energy functional is obtained by 
taking a trial state p® N in the TV-body free energy, which 
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yields 



B. Properties of the mean field density 



£ mf [p] = [ 



P W N> ; 



2 11 p{z)\og\z-z'W) 

/r 2 xr 2 

+ N- 1 [ p\ogp (VI.5) 



with 



W N , m (z) = \z\ 2 ~2-\og\z\ 



It has a unique minimizer g MF among probability mea- 
sures on R 2 . The basic result concerning its relation to 
the scaled single-particle probability density 

P ( n( z )= p N (z,z 2 , . . . ,z N )dz 2 ■ ■ .dz N (VI.6) 

JR2(N-D 

is as follows: 

Theorem 4 (Plasma analogy for QH trial states) 

There exists a constant C > such that for large enough 
N and any smooth function V on M 2 

($(z)- g UF (z))v(z)dz 

< ^(logJV/JVjVailVVlliacR.) + CiV- 1 / 2 ||Vy|| L oo (R2) 

(VI. 7) 

if to < A^ 2 , and 

^ } (z)-g MF (z))v(z)dz 

< C7V- 1 / 2 m- 1 / 4 ||y|| i ^ (R2) (VI.8) 

ifm > N 2 . 



The different form of the estimates (|VI.7j) and (|VI.8j) 
are due to different methods of proofs in the two differ- 
ent regimes and already hint at the fact that something 
special happens in the regime m oc N 2 , as we discuss in 
more details in Section fVl Bl below. 

Essential ingredients in the proof of Theorem U are 

• 2D versions of two classical electrostatic results: 
Onsager's lemma, and an estimate of the change 
in electrostatic energy when charges are smeared 
out (see e.g. [H|, Ch. 6). 

• the variational equation associated to (jVI.5l) and 
the positivity of relative entropies. 

These put together give precise estimates on the JV-body 
free energy functional. Using also the Cszizar-Kullback- 
Pinsker inequality and a new refinement of Onsager's 
lemma one eventually deduces estimates on the one-body 
density. We refer to 2l| for more details. 

In order to apply the theorem to estimate the energy 
in a potential like V Ut k{z) — W M 2 + k\z\ 4 , that increases 
at infinity it is necessary to study also the decrease of 
the mean field density at infinity to compensate for the 
growth of Vu t k- 



Dropping the last term in (|VI.5I) (the entropy term) 
one obtains a simplified functional 



pw Ni , 



p(z)\og\z-z'\p(z') (VI.9) 



that can be minimized explicitly. In fact, using the 
variational equation of the minimization problem and 
electrostatic arguments, in particular Newton's theorem, 
it is not difficult to prove that the minimizer, denoted 
by @ , takes the constant value l/(27r) in an annu- 
lus with inner radius R~ = W m/N and outer radius 
R+ = y/2+ (m/N), and is zero outside the annulus. In 
particular, for m — the support is a disc with radius 
\[2. Note that the length scales here refer to scaled vari- 
ables, cf. (|VI.2|) . In the original variables the lengths 
have to be multiplied by \f~N and the density by TV -1 . 

The reason for the label "el" in (|VI.9|) is that the func- 
tional and the shape of the minimizer g cl is entirely de- 
termined by the electrostatic terms in the free energy 
functional. This minimizer can be shown to approximate 
g MF in the metric that is defined by the inverse of the 2D 
Laplacian: Define for real valued functions (or measures) 
/ and g 



D(f,< 



f(z)\og\z~z'\g(z')dzdz'. 



(VL10) 

Then D(f, f) > if / / = with equality only if / = 0. 
The approximation result for g MF in terms of g cl is 



D{Q MF - Q el ,Q MF - Q el ) < CN-\ 



(VI.11) 



An analogous estimate of the difference between and 
g MF lies behind (|VI.7|) . Combining these two estimates 
one sees that (|VI.7j) holds with g MF replaced by g cl . 

From the formulas for the inner and outer radii 
it is clear that the shape of the g cl depends essentially 
on the ratio m/N . As long asmCW the support is a 
disc of radius 0(1) with a small hole around the center. 
As the ratio increases the hole becomes larger and for 
to 3> N the support is an annulus of mean radius R ~ 
yJm/N > 1 and thickness R + - R~ ~ ^N/m < 1. The 
exponential decrease of g MF for large distances r from the 
origin sets in for r — R + ^> max{i? + — R~ , A^ -1 } (again 
we are stating the results in terms of scaled variables, cf. 
(IVL2l) V 

The length scale A^ _1 is the coefficient in the entropic 
term of the free energy functional and it becomes larger 
than the length scale R + — R~ when to > A^ 2 . In fact, 
under this condition it is better to approximate g MF by 
dropping the interaction term D(p,p) from the energy 
functional, defining 



£ th [p] = P W N . 



N~ 



plogp (VI. 12) 
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where 'th' stands for 'thermal'. The minimizer of this 
simplified functional, denoted g th , can be computed ex- 
plicitly and is given by the density of the one-particle 
state (7rm!)- 1 /V n exp(-|z| 2 /2), i.e., 

g th (z) = (^m!)- 1 |z| 2m exp(-|z| 2 ). (VI.13) 

In the radial variable r = \z\ this density is approximately 
Gaussian around r = y/m with width 1 / \fm. 

This transition from an essentially electrostatic to an 
essentially thermal behavior of the minimizer of the 
mean-field functional (|VI.5|) and thus, via (IVI.7|) . of the 
one-body density of Laughlin's quasi-hole when m ^> N 2 
does not seem to have been noticed before. Indeed, as far 
as we know, the usual approximations in the literature 
have been based on purely electrostatic arguments, that 
is, implicitly, on the simplified functional (IVI.9[) . Our 
analysis shows rigorously that this approximation is cor- 
rect as long asm< N 2 , but breaks down in the opposite 
case where one should rather consider (IVI.12I) . 

Physically, the reason for this is that to minimize 
(|VI.5|) . the density has to be squeezed as tightly as pos- 
sible around the potential well of T4V,m- There are two 
limiting effects for the squeezing, one associated with the 
electrostatic term in (|VI.5I) and the other with the en- 
tropy term. Each term comes with its own length scale 
and the true density is spread over the larger of the two, 
which is the 'electrostatic' length */ N/m for m <C N 2 
and the 'thermal' length for m ^> A 2 , indicating a tran- 
sition in the regime m oc TV 2 . 

Computing the potential energy with the density g cl 
and optimizing over m leads to the same optimal value 
m op t as in (|V.8|) . In particular, the Laughlin state, i.e., 
m = 0, leading to a full disc instead of an annulus, is 
favored for u> > —2kN. Moreover, together with the 
estimates on the decay of g MF away from the support of 
g el , the energy upper bound pV.10[) is improved to 



E < ujN 2 + -kN 3 
3 



(VI. 14) 



and pVTTl) to 

u 2 N 



En < 



u . +(|*JV 3 + |mJV)(1. + ( .(1II (VI.15) 



at least for \ui\/k < TV 7 / 5 In A. (This limitation is due 
to the fact that the the current estimates of the decay 
of g MF are less than optimal.) Corresponding improve- 
ments hold for the angular momentum estimates: Eq. 
(|IV.1|) is replaced by 



La < (4/3) 1/2 A 2 



(VI.16) 



and the V3 in Eqs. (|IV.2I) and (|IV.4|) can be replaced by 
1/V3. 

Similar improvements of the angular momentum and 
energy estimates can be provided in the thermal regime 
using the density g th , but again due to lack of precision 
of the decay estimates, they can only be proved with our 
method when \oj\/k >■ TV 10 / 3 . 



VII. DISCUSSION 

We have studied a rotating Bose gas in a quadratic plus 
quartic trap where the rotational frequency can exceed 
the frequency of the quadratic part of the trap. Through 
the analysis of trial states for energy upper bounds and 
simple lower bounds we have obtained criteria for the 
ground state to be strongly correlated in an asymptotic 
limit. The lower bounds, although not sharp, are of 
the same order of magnitude as the upper bounds. The 
optimal trial state changes from a pure Laughlin state, 
with essentially constant density in a disc, to a modi- 
fied Laughlin state with a 'hole' in the density around 
the center when u> is negative and |w|/2A: exceeds N. As 
\uj\/k increases, the density becomes concentrated in an 
annulus of much smaller width than its radius. Within 
the annulus the density profile is still approximately flat 
until another transition occurs for \ui\/k of the order N 2 
when the density profile becomes approximately Gaus- 
sian in the radial variable. 

We believe that these results are sharp in the sense 
that our trial states, the Laughlin quasi-holes (Laughlin 
times giant vortex), are good approximations to the true 
ground states of (|III.5|) in the fully correlated regimes 
described in Theorem [2j We have no proof of this for 
the moment: Rigorous lower bounds to the energy that 
would match our upper estimates would require a better 
understanding of the properties of general fully correlated 
states. 

Finally we discuss the connection between our ap- 
proach and the experimental proposal of (Tlj . There it 
is argued that a convenient way of creating the Laugh- 
lin state in a cold Bose gas would be via a dynamical 
procedure. As in our approach an anharmonic tra p is 
used, but this only in a first step. The main idea of [1J| 
is to first use a "Mexican hat" anharmonic potential to 
create a rotating BEC with a large angular momentum 
concentrated in a giant vortex at the center of the trap. 
In a second step one would then stop the rotation and 
change adiabatically the trap from Mexican hat like to 
purely harmonic. With a careful tuning of parameters, 
one then ends up with a gas where interactions dominate 
the physics and thus strongly correlated states should 
emerge. Since the angular momentum is conserved dur- 
ing the adiabatic evolution, if one manages to first create 
a condensate with momentum N(N— 1), the final evolved 
state should be the Laughlin state, and [ll[ mostly em- 
phasizes this case. If one instead starts from a BEC with 
a larger momentum N(N — 1) + mN it is likely that the 
final evolved state will be the corresponding Laughlin 
times giant vortex (IVI.1|) . 

Indeed, as explained in Section [II A[ an interpretation 
of the Laughlin state is that vortices are bound to par- 
ticles in order to decrease the interaction energy. In this 
perspective, one may interpret the adiabatic evolution 
proposed in 11| as a redistribution of the vortices, all 
located at the center of the trap in the initial BEC state. 
During the evolution, vortices are bound to particles un- 
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til the interaction energy is zero. If one starts with more 
vortices than necessary to achieve this, i.e. with an angu- 
lar momentum strictly larger than N(N— 1), the 'excess' 
vortices will probably have a tendency to stay at the cen- 
ter of the trap to conserve rotational invariance, which 
leads to the form (|VI.1|) as a guess for the final evolved 
state. A possible way to further enhance this effect and 
stabilize the giant vortex could be not to completely turn 
off the anharmonic trap, so that the final trap also has 
some Mexican hat shape that would help pin the excess 



vortices at the origin. 
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